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Abstract: The scattering theory of Lax and Phillips, designed primarily for hyperbolic 
systems, such as electromagnetic or acoustic waves, is described. This theory provides 
a realization of the theorem of Foias and Nagy; there is a subspace of the Hilbert space 
in which the unitary evolution of the system, restricted to this subspace, is realized as a 
semigroup. The embedding of the quantum theory into this structure, carried out by Flesia 
and Piron, is reviewed. We show how the density matrix for an effectively pure state can 
evolve to an effectively mixed state (decoherence) in this framework. Necessary conditions 
are given for the realization of the relation between the spectrum of the generator of the 
semigroup and the singularities of the S- matrix (in energy representation). It is shown 
that these conditions may be met in the Liouville space formulation of quantum evolution, 
and in the Hilbert space of relativistic quantum theory. 

1. INTRODUCTION 

The unstable quantum system is an important example of irreversible phenomena in 
nature. Such systems, ranging from excited atomic states to short-lived elementary parti- 
cles, are characterized by what is generally observed to be an irreversible evolution. These 
phenomena raise the question of explanation of such processes from first principles. More- 
over, since most of the decay processes are observed experimentally to obey an exponential 
decay law, one expects this behavior to follow from very general assumptions. 

Irreversible evolution in the quantum theory has been described by the addition of non- 
Hermitian terms to the Hamiltonian, such that it has complex eigenvalues, and the induced 
evolution is non-unitary. Structures of this type were originally introduced by Gamow [1] 
who studied the effect of assigning complex eigenvalues to the energy spectrum, and hence 
introduced a kind of generalized eigenvector. Wu and Yang [2] parameterized the K- 
meson decay in this way. In this method, the non-Hermitian terms in the Hamiltonian are 
introduced phenomenologically, and may only indirectly be associated with some known 
interaction in a more fundamental Hamiltonian. 

Weisskopf and Wigner [3], in a well known paper in 1930, introduced an alternative 
approach to the decay problem. According to their approach, the evolution takes place 
in a Hilbert space which is a direct sum of two subspaces: the subspace of the decay- 
ing states and that of decay products. These two subspaces are stable under the "free" 
evolution induced by Hq, but are combined linearly under the full evolution induced by 
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H = H + V. In this Hilbert space, the evolution is unitary, and hence its generator, i.e., 
the Hamiltonian, is self-adjoint. The decay is described as the probability flow from the 
subspace of the decaying states to its complement, the subspace of the decay products. 
They studied perturbatively, for the single-channel case, what has become known as the 
survival amplitude 

A(t) = (V,e- iH V), (1-1) 

which is the probability amplitude for the system to remain in the discrete state until time 
t. In the following we will describe this approach, and pose critical problems, motivating 
the development of a more general theory. 

Let us denote the projection operators on these two subspaces as P and P, such that 
P + P = 1. For the decay problem, the basic quantity is the reduced motion 

U'{t) = PU(t)P. (1.2) 

where U(t) = e~ lHt , which governs the time evolution of the subspace P7i of the unstable 
states. From this one can derive the decay law of the unstable states. If {<pi} is an 
orthonormal basis of PH, the probability that an unstable state <p, which exists at time 
t = 0, is in the subspace PH of unstable states at time t is given by 

p(t) = J2 U(t)<P)\ 2 = Tr(U'(tyU'(t)P^ (1.3) 

i 

where P^ = \4>){4>\- 

The total evolution operator U(t) = e~ xHt , and the resolvent R(z) = (z — H)~ x are 
related to each other by the (inverse) Laplace transform 

U(t) = ^- j) R(z)e~ lzt dz, (1.4) 

where the integration contour is around the spectrum of H. If we project this operator 
into the subspace PH., we can obtain a similar relation which expresses the reduced motion 
U'(t) in terms of the reduced resolvent R'(z) = PR(z)P: 

U'{t) = — I R'(z)e- lzt dz. (1.5) 
2ni J 

By differentiating Eq. (1.3) and setting t = 0, one sees that the initial decay rate 
is necessarily zero (providing that the Hamiltonian is defined on the initial state); in 
fact, it is easy to show that the change in p(t) is 0(t 2 ). The intermediate and long 
time behavior follow most simply by an examination of the relation (1.5). Deforming the 
contour of integration which runs below the real positive spectrum of H to the negative 
imaginary axis, where its contribution for large times is small, the remaining contribution 
of the contour running above the spectrum of H can be estimated by bringing this contour 
continuously through the cut. When resonances exist, this contour will pass through simple 
poles on the way to the negative imaginary axis in the second sheet (as can be explicitly 
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demonstrated [4] in the Lee-Friedrichs model [5]). The residues of the poles may dominate 
the time dependence for intermediate times, and give the approximate (due to the presence 
of residual contributions from the integrals along the negative imaginary axis) exponential 
decay behavior. For very long times, the pole contributions disappear, and the remaining 
integration around the branch cut results in an inverse power law asymptotic behavior [6] . 

It is not difficult to see that an irreversible evolution must be described by a semigroup 
[7] (for the reversible case this is a group induced by a unitary transformation), where we 
define a semigroup as follows: 

Let {Z(t)} be a family (over t) of operators on a Hilbert space; then Z(t) is an element 
of a semigroup if 

Z{t x )Z{t 2 ) = Z(h + 1 2 ), h, t 2 >0 (1.6) 

The semigroup is said to be strongly contractive if ||Z(t)|| — > 0, for t — > oo, where \\A\\ is 
the operator norm of A. On the other hand, it can be shown that the reduced motion, as 
described above, cannot generate a semigroup [8]. 

There is, furthermore, another, perhaps more fundamental problem associated with 
the general method of Wigner and Weisskopf; this is that the expression (1.1) for the 
survival amplitude implicitly assumes the existence of a linear superposition (we restrict 
our discussion here to the one-channel case) 

e-^V = A(t)iP + X (t) , (1.7) 

where represents the decayed system and (ip,x(t)) = 0. In general this linear su- 
perposition does not correspond to any physical situation in our experience; a short-lived 
particle, for example, is seen as either the particle before the decay, or the decay prod- 
ucts at a certain time, which can not be predicted. This linear superposition does not 
correspond to the object that we see experimentally in such a process. 

In the framework of the theory of Weisskopf and Wigner, techniques have been devel- 
oped which are capable of displaying the exact semigroup behavior of an unstable system 
[9] . As we have remarked above, the Lee-Friedrichs model [5] provides a simple but useful 
example of an unstable system for which the evolution equations are completely soluble 
[4]. The eigenvalue equation Hf(z) = zf(z), for / e Tl has a formal solution which does 
not, however, satisfy the equation. The condition to satisfy the equation coincides with 
the condition for a pole in the resolvent R'(z), and can only be satisfied in the second 
Riemann sheet. The scalar product of the eigenvalue equation with a vector for which the 
unperturbed energy representation is analytic in a domain containing the second sheet pole 
can be continued to the second sheet, and at the pole position, the equation is identically 
satisfied. Since the "eigenvector" obtained in this way is in the space dual to a set of 
vectors with this restrictive analyticity requirement, corresponding to a subspace of H, it 
is an element in the large space of a Gel'fand triple (rigged Hilbert space). The definition 
of this vector depends on the domain of analyticity chosen, and its physical interpretation 
is not clear, except for the fact that its (extended) unitary evolution is that of an exact 
exponential decay. 

The problems discussed above are essentially related to the attempt to describe an 
unstable system in a framework more suitable to the description of reversible phenomena. 
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In what follows we will show another approach to irreversible phenomena which attempts 
to solve these difficulties. 

2. LAX-PHILLIPS THEORY AND THE EXACT SEMIGROUP 

The characterization of a system undergoing an irreversible process cannot, in prin- 
ciple, be specified at a given instant of time. In fact, the physical quantities describing 
such processes involve time measurements (that is, measurements of the time at which 
certain defined phenomena occur). Therefore, the information about the decay which is to 
be deduced from the state is associated with its distribution in time, an essential property 
of the system, just as the location or momentum of a quantal particle. The time variable 
is, from this point of view, an internal degree of freedom of the system, which provides a 
framework for the description of interactions which can influence the structure of the state 
[10]. The dynamical evolution of the system involves a change in its internal structure, 
including its distribution in t along with other observables characterizing the state. This 
evolution, parameterized by the laboratory time r (which is not a dynamical variable), is 
defined on a Hilbert space H with t in its measure space with norm given by (e.g., with 
Lebesgue measure) 

J uifdt = (2.1) 

where the norm in the integral is taken as the norm in 7it, a member of a family of auxiliary 
Hilbert spaces (all isomorphic), defined for each t. 

The theory of Lax and Phillips [11], designed for systems of hyperbolic differential 
equations describing the scattering of, e.g., electromagnetic or acoustic waves, and the 
Floquet theory [12] for periodic time dependent quantum mechanical problems are exam- 
ples of such a structure. Piron [7] has shown that methods of this type are applicable 
to the general time dependent quantum mechanical problem. Recently, Flesia and Piron 
[13] have shown that scattering problems in quantum theory can be put in the form of 
Lax-Phillips theory (Horwitz and Piron [14] have discussed its applicability to the problem 
of the unstable system) by forming a direct integral of the quantum mechanical Hilbert 
spaces Ht over t in order to construct a larger space H which includes t in its measure 
space. 

Lax-Phillips theory [11] assumes the existence of a one-parameter unitary group of 
evolution on a Hilbert space H, and incoming and outgoing subspaces V- and V + such 
that 

U(t)V+ C V+ , for all r > 
U(t)V_ C £>_ , for all r < 
f]U(r)V ± },= {0} 



\JU(t)V± = H (2.2) 

r 

where r is the evolution parameter identified with the laboratory time. It follows from a 
theorem of Sinai [15] that H can be foliated in such a way that it can be represented as a 
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family of (auxiliary) Hilbert spaces in the form L 2 (— oo, +00; Tit), over Lebesgue measure 
in t , and all the Tit are isomorphic (we therefore sometimes refer to these spaces simply 
as 7i) and determined up to unitary equivalence. The scalar product in H is given by 

/oo 
(ft,9t)n t dt. (2.3) 
-co 

Lax and Phillips show that there are unitary operators WZ l ,WZ X which map the 
elements of H into representations, called the outgoing and incoming translation represen- 
tations, for which the evolution is translation in t. The subspaces T> + ,T>_ correspond to 
the sets of functions with, in these representations, support in semi-infinite segments of 
the positive and negative t-axis respectively. They define the S matrix abstractly as the 
map from the incoming translation representation to the outgoing one, i.e., £ = W^VF-. 
This map is defined up to unitary transformations on the auxiliary spaces {Ht}, and refers 
to the equivalence classes for which the incoming and outgoing representations have the 
property that the evolution is represented by translation. 

Lax and Phillips furthermore define the operator 

Z(r) = P+U(r)P. (2.4) 

on 7i, where P± is the projection on the orthogonal complement of T>±. This operator 
vanishes on V± and maps the subspace 

K = HQ (V+®V_), (2.5) 

into itself. These mappings form a semigroup [11], i.e., for t\, r 2 > 0, 

z(n)z(T 2 ) = z( Ti + t 2 ), (2.6) 

and this semigroup is strongly contractive, i.e., for each 4> e /C and any e, there exists a 
such that 

\\Z{r)4>\\n < c (2-7) 

for t > Tfj). It can be shown that Z{t) is just the unitary evolution U(t) projected into 
the subspace /C. Since the states which lie in the subspaces T>±, in the case of scattering, 
describe the incoming and outgoing waves which are not influenced by the interaction, 
the states which lie in K, describe the unstable states, i.e., resonances of the scattering. 
From this point of view, the Lax-Phillips semigroup is analogous to the reduced motion 
discussed in the previous section. 

This theory constitutes a constructive realization of the theorem of Foias and Nagy 
[16], which states that given a semigroup on a Hilbert space H, there is a bigger Hilbert 
space H which contains it, in which the evolution is a one-parameter unitary group, and 
this unitary group restricted to H is that semigroup.* 

* The theorem states that this construction is minimal. We conjecture but have not 
proved that the Lax-Phillips construction is minimal. We thank G. Emch for a discussion 
of this point. 
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Flesia and Piron [13] have shown that the quantum theory may be embedded in a Lax- 
Phillips theory by considering the family of Hilbert spaces of the usual quantum theory 
on the parameter t as the auxiliary spaces of Lax and Phillips; the large Hilbert space 
H is then the direct integral of these quantum mechanical spaces over all values of the 
time t with Lebesgue measure. The form of the theory adopted by Flesia and Piron [13] 
distinguishes the elements of these equivalence classes, and constructs an S-matrix which 
maps the auxiliary space in the incoming translation representation to the auxiliary space 
of the outgoing one. In the model that they use to illustrate this structure, this map 
corresponds to a pre-asymptotic form of the S-matrix of the usual scattering theory. Their 
model assumes that the subspaces X>+, £>_ are represented in the "free" representation, for 
which the free evolution is translation, by L 2 (— oo, p_; 7i), L 2 (p + , oo; Ti), respectively. In 
the limit in which the interval between the two semi-infinite regions of support tends to 
infinity, their S-matrix becomes the usual S-matrix. In this construction, Flesia and Piron 
assume the form 

VtV = W t (r)^, (2.8) 
where, since Wt(r) represents an evolution, it follows that 

Wt+^Wtin) = Wtin + r 2 ) . (2.9) 

Lax and Phillips prove that the S-matrix (in their construction) is a multiplicative 
operator in the spectral representation of the generator of the unitary evolution K (which 
is the Fourier transform of the translation representation), i.e., 

{Sip) a = S(cx)Vv, 

and that the eigenvalues of the generator of the semigroup Z(t) correspond to the singu- 
larities of the analytic continuation of S(cr). The eigenstates corresponding to these eigen- 
values are analogous to the generalized eigenstates found in the framework of Weisskopf 
and Wigner, as discussed in Section 1. Thus, the S'-matrix contains all the information 
about the unstable states . It can be seen [10], however, that the S matrix obtained from 
a model in which the evolution is given in the form (2.8) has no t-dependence, and hence 
its spectral representation is trivial. In this form, one therefore has no relation between 
the singularities of the S'-matrix and the spectrum of the generator of the semigroup. 

Although the generalization of Lax-Phillips theory by Flesia and Piron [13] provides a 
new point of view for scattering theory, we see that to extend the theory further to include 
a description of the evolution of an unstable system, it is necessary to generalize the law 
of evolution to that of a nontrivial integral operator over the time. 

The most general linear evolution law has the form 

/+oo 
W t , tl {r)^dt' . (2.10) 
-oo 

We shall show that this type of evolution, which goes beyond the formulation of Flesia 
and Piron [13] and Floquet theory [12], can correspond to unitary evolution in H with a 
nontrivial S-matrix for which the singularities of its Fourier transform are associated with 
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the spectrum of the generator of the Lax-Phillips semigroup. As we shall show below, the 
form of the evolution law (2.10) has a natural realization in Liouville space as well as in 
the framework of relativistic quantum theory. 

Let us now study for this general evolution, some properties of the S-matrix 



(Sty)t = J S t ,t^t>dt', 



and show that in this general case the S-matrix must have the form St,? = S(t — t'). 
Using the definition S = W_jT W_, where 



W± = s- lim U(-t)U (t), 

T— >±0O 



we find 



But, 



(SiP) t = s- lim (^o(-ri)^(n)?7(r 2 )?7o(-r 2 )^) t 

Tl,T2—HX> 

(C/o(-ri)C/(T 1 +T 2 )t/ (-T 2 )V;) t = {U{ Tl + T 2 )U Q {-T 2 )^)t +Tl = 
= |f 1+Tl ,t'(n+r 2 )(f/o(-r^Wi' = 

= J Wt+ T1 ,?(T 1 +T 2 )lp?+ T2 dt' = J Wt+ Tl ,?- T2 (T 1 +T 2 )lp?dt', 

and therefore the matrix elements of S are 

S t ,? = s- lim W t+Tl ,?-T 2 {Ti + T 2 ) = 

T\,T 2 — >00 

= s - # lim W t -t>+Ti,-T>(T[ + t' 2 ) = S(t - (2.11) 



(where t[ = T\ + t' t' 2 = t 2 — t' ) . This is a very important property of the S'-matrix, 
according to which, when one goes to the spectral representation ifj a = f e~ iat ifj t dt, the 
S'-matrix takes the simple form 

S(j,(j f = ^ J e-^St^e^'t'dtdt' = 6(a - <j')S(a) 

where 

S(a) = J e- iat S{t)dt, (2.12) 

i.e., in this basis the S-matrix is diagonal, and the S-operator is multiplication on the 
subspaces, labeled by a, of {H a }, the set of (isomorphic) Hilbert spaces which are the 
Fourier dual to the set {H t }. This result can be obtained also by looking at the definition 
of the S-matrix, 

S = s - lim Uoi-nMn + t 2 )u {-t 2 ) 

Ti,T2— >oo 
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from which it follows that 

SU (r) = U (t)S. 

Since Uq(t) is the translation operator one obtains the result [S , id t ] = (which corre- 
spond to the usual result of scattering theory [ S , H ] = 0) . It follows from this commu- 
tation relation that S t ,t' = S(t — t'). 

We show now that under the general evolution (2.10), the semigroup is contractive. 
Let us calculate the generator of the semigroup B of Z{f) = P + U(t)P_. We use the 
free translation representation in which both T>± have definite support properties. In this 
representation, 

Z(r) = P+U(r)P. = E(p)U(r)(I-E(0)), (2.13) 

where E(t) is the spectral resolution corresponding to Tq, the free-time-operator (the 
conjugate of K which is, in the free translation representation, — id t ). Then, the generator 
(in the subspace K) of Z{f) is 

B — i lim ^ ~ fc = i lim = ^ = g <°» = fc = 

= E(p)K(I - E(0)) = P + KP_. (2.14) 

According to the requirements on T>± the matrix elements of k, the self-adjoint kernel over 
t, t' contained in K distinct from the t-derivative [10], between states from D_ to T> + , or 
T>± to /C vanish, and therefore 

B = P+K P- + k k ■ (2.15) 
An operator B is called dissipative [17] [18] if 

<0, (2.16) 

for all <p i n the domain of B. Since kx, is self-adjoint only the first term determines whether 
the operator is dissipative, i.e., this property does not depend on the perturbation. As 
shown by Horwitz and Piron[14], the operator P + KqP_ is, in fact, dissipative. It is known 
[18] that Z(t) is a contractive semigroup if and only if its generator is dissipative. It 
therefore follows, independently of (self-adjoint) interaction, that the semigroup Z(t) is 
contractive. We see from this [14] the essential mechanism of Lax-Phillips theory. The 
non-self-adjointness of P+KqP- corresponds to the restriction of —idt to a finite interval, 
so that, in fact the operator has imaginary eigenvalues. In the presence of interaction 
(non-trivial k) : these eigenvalues emerge as the actual eigenvalues of B, corresponding to 
the singularities of S(a). 

We remark that the direct integral space provides a framework as a functional space 
for quantum mechanics in which the Nagy-Foias construction can be realized, i.e., for which 
unitary evolution can be restricted to a contractive semigroup. We shall now introduce an 
extension of the conceptual framework which considers the set {?/>t}, corresponding to the 
Lax-Phillips vector ip, as an ensemble of the same type, for example, as {ip(x)} G H, where 
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x is a point of the spectrum of the position observable, in the usual form of the quantum 
theory. In concluding this section, we investigate some consequences of this interpretation. 

In particular, we discuss some properties of the time operator and the realization of the 
superselection rule in time. In the next section, we discuss the possibility of decoherence 
in H induced by the unitary evolution in H. 

There are three distinct types of time operator. One, which we call the incoming time 
operator T m , provides a spectral family in terms of which the incoming representation can 
be constructed, and in which functions in T>- have definite support and functions in Ti 
evolve by translation. In this representation, the norm of the evolving states in (— oo, 0) 
must decrease. After sufficient laboratory time r passes, the states evolve to V +: and in 
the outgoing representation, provided by the spectral family of the outgoing time operator 
rpout^ thgy h ave definite support in (p, oo). The mapping of functions in the incoming 
representation to the outgoing representation is provided by the Lax-Phillips S'-matrix, 
and the time operators are related by 



The third type of time operator corresponds to the "free" representation and is related to 
rpmrpout ^ e L ax _phillips wave operators. The spectral family for this operator provides 
the "standard" representation (analogous to Dirac's choice of "standard" spectral families), 
which we have used above. 

There is an interval, in general, when the system is in interaction, and its state is 
neither in £>_ nor T> + . The expectation value of the operator T in in the state t(j T projected 
into /C ®T> + (corresponding to the projection P_) can be interpreted as the interaction 
interval. If the system in interaction is considered as an unstable particle (a resonance), 
this interval is its age after creation at t = 0. The expectation value of T tn then moves 
out of (—00, 0). The expectation value of T in in the state P_ip T is 



here, | i n (t\P- r ip T ) | 2 is the probability density for the age t at time r, an intrinsic dynamical 
property of the system. The positive value that the expectation value develops corresponds 
to the average age. One can similarly compute the expected time after decay, the expected 
lifetime, and the expected value of any other observable of interest as a property of the 
unstable system. 

We then understand the subspace /C as corresponding to the unstable system. 

The structure of the theory is somewhat similar to the Wigner-Weisskopf idea, in that 
a subspace is associated with the decaying system. The decay of the system is associ- 
ated with the probability flow out of the subspace. As in the original Wigner-Weisskopf 
formulation, the process of decay may be represented as a continuous evolution from the 
original unstable state to the final state through a changing linear superposition. In this 
framework, let us choose a vector tp in the subspace /C to represent the state of an unstable 
system. Then, under the full evolution, 



T 



tout 



ST in S \ _ 



(2.17) 




(2.18) 



= (V>,2(r)V), 



(2.19) 
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so that the reduced evolution is an exact semigroup. 

Moreover, in the Lax-Phillips theory the expectation value of an observable which 
is decomposable in the free or outgoing representations, where T> + has definite support 
properties, necessarily reduces to the sum of the expectation values in the subspaces K@V- 
and in the subspace V + (the decay products), i.e., 

(A) = f aw*, am = e / dt ^> ( 2 - 2 °) 

•* M=D ± ,K 

Note that there are no cross terms. There is, therefore, an exact supers election rule 
for measurements of the system by means of such decomposable operators. 

3. APPLICATIONS 

a. Measurement according to Namiki and Machida 

Recently, Machida and Namiki [19] have proposed a measurement theory based on a 
direct integral space of continuously many Hilbert spaces and a continuous superselection 
rule. As pointed out by Tasaki et al [20] , although they had some success, their theory has 
a conceptual difficulty. Indeed, in their theory, while the apparatus is described by many 
Hilbert spaces, the system corresponds to a single Hilbert space as in the conventional 
theory. Thus, one needs to specify the boundary between the system and the apparatus. 
As discussed by von Neumann, this is impossible. 

Most measurement processes are concerned with measurements of observables which 
are time-independent in the Schrodinger picture. Therefore, if two different Lax-Phillips 
states give the same expectation value for all time-independent observables, these two 
states are essentially indistinguishable. In this sense, we define the following: 

1. A Lax-Phillips vector ip G Tt is called " effectively pure" if there exists a pure state 

Po = 4>o4>o, 0o g H, 

such that 

(i)v, =Trp A= (0o,^40o), (3.1) 

where A is the "lift" of A on 7i to 7Y, for every element of the algebra of bounded linear 
operators associated with the spectral families of the time-independent observables* 
on the original space T~L. 

2. A Lax-Phillips vector is called "effectively mixed" if no such (pure) po exists. 

It can be shown[10] [20] that ip = {ip t } G H is effectively pure if an only if it has the 

form 

i>t = f(t)<fo. (3.2) 



* We wish to emphasize that what is meant is explicit time-dependence in the 
Schrodinger picture; we do not refer here to the dynamical time-dependence that may 
arise in the Heisenberg picture if A is not a constant of the motion. 
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We now discuss the possibility of decoherence, or the evolution from effectively pure 
to effectively mixed states. First, we consider the Schrodinger evolution for a time- 
dependent Hamiltonian. The solution of the time-dependent Schrodinger equation can 
always be written formally as tpt = U(t,t')ipt', where U(t,t') satisfies the chain property 
U(t, t')U(t', t") = U(t, t"), and can be expressed in terms of the integral of a time-ordered 
product. We define W t (r) = U(t + r, £), and lift the evolution to H as follows 

Vt + r = Wt(T)Vt, (3-2) 

where W t (r) is given by (T implies the time-ordered product) 

Wt(r) = T (e~ l $t +T H{t ' )dt '^j . (3.3) 
For this kind of time-evolution we obtain 

(i)^ = Jdt (W t (T)ip t ,AW t (r)ip t )n, (3.4) 
where we have taken the normalization as unity. For the effectively pure states we have 

(A)* = Jdt If^WWtir^AWtir^n. (3.5) 

It follows from our previous argument that the effective state corresponding to (3.5) is 
mixed-like if W t {r)4>Q ^ W t '{r)4>Q (i.e., the state p^ induced from ipJ +T = Wt(r)ipt = 
f(t)W t (T)(j) is not pure in H). This result is true for the generalized evolution W U ' of 
(2.10) as well. 

b. Intrinsic decoherence in classical and quantum Liouville evolution. 

It has long been emphasized by Prigogine and his co-workers [21] that the natural 
description for the evolution of a system with many degrees of freedom is that of the 
evolution of the density matrix p, through the Liouville equation, 

i^ = {H,p}. (3.6) 

The density matrix p (p > 0, Trp = 1) has the property that Trp 2 < 1, where the equality 
is attained only for a pure state. In general, one considers the space of Hilbert- Schmidt 
operators A for which 

TrA*A < oo; (3.7) 

the positive (normalized) elements of such a space correspond to the physical states, the 
density matrices. On this space, the commutator with the Hamiltonian H defines a linear 
operator C , called the Liouvillian, for which 

i% = £p, (3-8) 
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where one assumes that C is self-adjoint in the Liouville space. 

The Hamiltonian evolution of states in classical mechanics is known by the Liouville 
theorem to be non-mixing, i.e., to preserve the entropy of the system [22]. The same 
property holds for the quantum evolution as well, and follows from the unitarity of the 
evolution operator. This has been an obstacle to the consistent description of irreversible 
processes from first principles [23]. The usual use of techniques of coarse graining or 
truncation to achieve a realization of the second law does not follow from basic dynamical 
laws, and is fundamentally not consistent with the underlying Hamiltonian structure [24] . 
We shall now show that the existence of a time operator in the Liouville space provides 
a natural and consistent mechanism for the decoherence of physical states, i.e., that pure 
states become mixed during the evolution, both for quantum and classical systems. 

In particular, for a Hamiltonian of the form of the sum of an unperturbed operator 
Hq and a perturbation V, i.e., H = Hq + V , the corresponding Liouvillian is 

£ = £ + £ I . (3.9) 

Now suppose we consider the "time operator" To, conjugate to Co (with spectrum 
(—00,00); it satisfies 

[T , £ ] = i. 
Then, in the spectral representation of To, 

o(t\[To,Co]\t'} = iS(t-t'), 

or 

(£-0 O (*|£o|0o = iS{t-t'). (3.10) 

It follows that 

o(t\C \t')o = -id t 5(t-t'). (3.11) 

Hence, 

o(£|£|£'>o = -id t 8(t-t') + (*|£/|£')o, (3.12) 

where the last term is, in general, not diagonal. 

The method that we have described above applies as well to the formulation of classical 
mechanics on a Hilbert space defined on the manifold of phase space which was introduced 
by Koopman [25] and used extensively in statistical mechanics [24] . Misra [26] has shown 
that dynamical systems which admit a Lyapunov operator necessarily have absolutely con- 
tinuous spectrum in (—00, 00); therefore one can construct a time operator on the classical 
Liouville space for such systems. The expectation value of a t-independent operator defines 
a reduced density function in the form 

J dt p t {p), 

where (3 is the set of variables remaining after extracting t as a function on the manifold of 
the measure (phase) space. Since a pure state is defined by a density function concentrated 
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at a point of the phase space, a state which is effectively pure must have the form 8{j3 — (3q). 
The equivalence class associated with this reduced density contains mixed states as well, 
such as p(t, (3) = S((3 — (3o)f(t) corresponding to a non-localized function on the phase 
space. The structure of the theory, and the conclusions we have reached, are therefore 
identical to those of the quantum case. 

c. Relativistic quantum mechanics. 

The form of relativistic quantum mechanics introduced by Stueckelberg [27], extended 
to the many-body case by Horwitz and Piron [28] , covariantly describes the evolution of a 
system according to the Stueckelberg-Schrodinger equation 

where M is an intrinsic property of the particle ("on-shell" mass). The classical form of 
this theory has for its Hamilton equations 

dx^_dKo_p^ 

dr d PfI M' 1 } 

and, therefore, eliminating dr, one obtains the standard relativistic relation 

Since the d'Alembertian, corresponding to the operator Kq, has spectrum (— oo, oo), there 
exists an operator £ which satisfies 

[K , £] = i. (3.16) 

Note that the operator t of the relativistic theory will not serve this purpose, since its 
commutator with Kq is iE/M, which only approaches i in the non-relativistic limit. 

If £ is a function of x, t, we may construct the transformation function (£', (3\x) using 
the defining commutation relation, i.e., 

(Z',P\K Z-ZK \x) = i{Z',P\x), 

or 

so that we obtain the defining equation [29] 

d d^d 

We thus see that the relativistic quantum theory provides a natural framework for 
the Lax-Phillips formulation of the description of an unstable system. It is interesting 
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that the continuous spectrum of K is essential to the construction; this implies that we 
must have both positive and negative mass-squared states in the spectrum, i.e., that the 
so-called tachyons, at least in the form of intermediate states, play a fundamental role in 
the relativistic description of unstable systems. 

One might ask how such a mechanism could survive in the non-relativistic limit. It is 
interesting to study this limit; even though the Galilean world is an idealization which is 
not realized physically, the velocity of light c is very large. To study this limit, we consider 
the condition 

E - Mc 2 = e < co, (3.18) 
for c — > oo, used, for example, in ref. [30]. We then define the variable m such that 

E = c^/p 2 + m 2 c 2 , (3.19) 

so that 



Defining 
we see that 



E-Mc 2 = (m-M)c 2 + mc 2 {Jl + ^ - lj. (3.20) 

V m 2 c 2 

7] = (m- M)c 2 ^ 0, 

2 

E = Mc 2 + V + ^ + f( V ,p 2 ), (3.21) 

an integral kernal on the wave functions, where the integral operator is 0(l/c 2 ) . 

The general structure of the relativistic Lax-Phillips theory therefore remains in the 
Galilean limit (for finite but large c) . The experimental signature of such an additional term 
in the Hamiltonian would be, for example, an interference effect in time. The experiment 
would be of the same design as the test of such an effect in the full covariant relativistic 
theory [31]. 
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